Merger Dynamics in Three- Agent Games 
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We present the effect of mergers in the dynamics of the three-agent model studied by Ben-Nairn, 
Kahng and Kim and by Rador and Mungan. Mergers are possible in three-agent games because 
two agents can combine forces against the third player and thus increase their probability to win a 
competition. We implement mergers in this three-agent model via resolving merger and no-merger 
units of competition in terms of a two-agent unit. This way one needs only a single parameter which 
we have called the competitiveness parameter. We have presented an analytical solution in the fully 
competitive limit. In this limit the score distribution of agents is stratified and self-similar. 
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I. INTRODUCTION 



The use of methods inspired from physical principles has recently been of wide utility in studying various systems. 
Models with explanatory and predictive powers have been applied to biological, social, political, economical, animal 
behavioral and many other phenomena. The common property of all these models is that they all involve a collection 
of agents interacting with a well defined set of rules. The set of rules may describe a competition for earning a game, 
for exchanging an attribute say money, energy or points, or to fight for a position in space. 

Recently such a model based on two-agent units of competition is introduced and was applied to sports dataQ-Q. 
Later that model is generalized to three-agent units providing qualitative understanding of emerging social structures 
That model being deterministic did cover the full range of possibilities. This extension is achieved in [f| where 
the full phase space of three-agent games is found yielding new social structures. 

In the present paper we study an intriguing extension of the three-agent model • the possibility that two agents 
could combine forces against the third one, in other words merger dynamics. In [j] and @ it was shown that the 
competitive subspace of three-agent games yields a continuous point gain rate for agents. Implementing mergers into 
the model we find that the agents condense at particular values of income rate. That is, as a result of mergers the 
society of agents becomes stratified. In the limit where the model becomes the most competitive the distribution 
for income rate becomes self-similar. 

We first review the main results of [5[ to introduce the three-agent model. The remaining parts of the manuscript 
is devoted to the study of merger dynamics. 



II. REVIEW OF THREE AGENT GAMES 



In this chapter we will review the dynamics of three agent games to have a structured manuscript. We shall omit 
various details as those were already discussed in depth in [5[. We shall however add some aspects that are not 
discussed in the mentioned paper and we shall add emphasis on points relevant for the next chapter where we will 
implement mergers into the model. 

Let us consider a collection of N agents with a given distribution of points. We pick three of them randomly and 
order their points say as L > M > S. We shall give one point to one and only one agent according to the following 
rules, 

• {L > M > S} =► {P, T, Q} 

{L = M >S}^{^f, £±£ Q} 



o 



o 



{L > M = S} => {P, ?±Q, T±Q] 



Here the lists on the right represent the winning probabilities of the teams with points listed on the left. In view of 
later chapters we refer to this model as the single rule model emphasizing the fact that there are no conditions on 
the points other than their ordering. As can be inspected, cases when some agents have equal points are resolved 
on the basis of equal likelyhood. That after every game one agent surely advances requires the normalization of the 
probabilities 



P + T + Q = 1 . 

Now let us denote the probability to pick an agent with point x at a particular time as f x , this should really be 
taken as a shorthand notation for f(x,t). After a competition some teams might leave this region towards x + 1 
and some teams might enter it from x — 1 by winning a competition in either case. This suggests the following local 
conservation law 



■£ = E f*-ifvfv' w ( x - x > w, y') - E f'fvfv w ( x > y> y') • w 

v,v' y,y' 

Here W(x, y, y') denotes the probability that the agent with point x will win againts two others with points y and y ' . 
The microscopic rules above completely define what W is. Furthermore, the right hand side is cubic in /. This is so 
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because we are picking three agents out of the collection and the probability to pick a team with a given point x is 

fx- 

Since 

E/» = 1 ' 

X 

it is immediate that (fTJ) also implies, as it should, the global conservation of the total number of teams, as can be 
checked by performing a sum over x. 

The time variable in ([!} has an arbitrary scaling which can be compensated by an overall factor in the definition 
of W since the former represents in essence a rate. The natural scale is such that the average points of teams is given 
by, 

°° t 

5(t) = 5> = 3 ' ( 2 ) 

x=0 

meaning that (on average) each team participates in a single game while we increment the time variable by one unit. 
As only one of the participating teams in a game wins and advances its score by one, equation ([2} follows. One can 
equivalently say that the average speed with which agents increase their points is 1/3. This normalization also means 
that the maximun theoretical point an agent could have acquired at time t is simply t. 

The presence of sums over the discrete indices on the right hand side of ([T]), results in a coupled set of differential 
equations. For the model at hand these can be much simplified by defining 



y=x* 



Here x* represents the smallest point below which there are no agents. Since agents do not loose points this value 
does not change in time and is defined by the initial point distribution as f x = at t = for x < x* . Without loss of 
generality on can take x* — and confine the points to positive values. Also from the definition of F x the following 
is immediate 



Summing (fTJ) over x we obtain 



dF x 
dt 



fx = F x+1 - F x . (4) 



-f x -iY,W(x-l,y,y')f y f y ,. (5) 



We note that 



fx-iY, w ( x -^y^') h fv (6) 

yields the probability that a team with score x — 1 will win any possible choice of single competition with two other 
teams. Working out the sum using the rules we get the final form of the master equation, 

dF 

=~U-i [pF^ + 9 (1 - F x ) 2 + 2iF x _i(l - F x )] 
o(P + *) f a p 

-Jx-l^x-l 
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Vl-x- (7) 
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The terms on the first line of ([7]) represent the bulk of interactions between three players with different scores. The 
second and third lines represent the cases where two players have identical scores x and the last term represents the 
case when all the teams have equal score. The effect of terms denoting units of competition where one picks some 
agents with equal score are irrelevant for the late time dynamics of the system as they die out in time. We will call 
these type of terms the interface terms and the rest as the bulk terms. So as time goes by, in a thermodynamic limit 
where the number of teams ranges to infinity, the majority of the contributions to the dynamics will be governed by 
the bulk terms. On the other hand as time goes by, almost every team will accumulate a certain number of points 
which, in general, will be larger than a single point. These considerations allow one to go to a continuum limit where 
the terms like F x -x are expanded in terms of the derivatives. A first order approximation, where one considers only 
the bulk terms, results in the following, 



dF 



.21 G'(F) = dG{F) 
dx dx 



(8) 



with 



G'(F) ee PF 2 + 2TF(1 — F) + Q(l - F) 2 . (9) 

We will refer to this approximation as the hydrodynamical limit. We would like to mention that with this limit we 
still have x = t/3, which corraborates the neglection of the interface terms. 

At this point we would like emphasize that the hydrodynamical limit for any model involving the competition of 
three agents will take the form 



^ = -^J J dydy ^w w[x ^ yh 

for a given function W(x, y, y'). 

In Q the solutions to (|5]) was studied via the method of characteristics, for all cases satisfying P + T + Q = 1, 
with the initial condition F(x, 0) = O(x) meaning that all agents starts with zero points. As is well known this is the 
Riemann problem. The asymptotic solutions was found to have the form F(x,t) — F(x/t) as it is usually the case 
with hyperbolic equations coming from conservation laws. In our normalization of the time variable the maximum 
possible theoretical point an agent can have at time t is simply t, thus z = x/t satisfies < z < 1. Furthermore one 
can show that x = t/3 yields 

/ dzF(z) = 2/3 . 
Jo 

In terms of the variable z — x/t one can recast @ into 

dF 

— [-z + PF 2 + 2TF(1 — F) + Q(l - F) 2 } = , (10) 
subject to the requirements, which we have discussed, that F has to satisfy. The solutions fall into the following 



categories 




Regime 


C~ : 


: for P > T > Q and T < 1/3, 


Regime 


C° 


: for T= 1/3 > q, 


Regime 


C+: 


: for T > 1/3 and P > T, 


Regime 




: for P < T and Q < 1/3, 


Regime 


s 


: for Q > 1/3 > P, 


Regime 


Cs: 


: for Q > T and P > 1/3. 
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FIG. 1. Representative solutions of the single rule model relevant to our discussion. For P = 1 we have an example of C~ solution 
which in this instance is simply F(z) = \fz. For T = 1 we have an example of Cg regime; the curve is F(z) = (1 — s/l — 2z) /2 
upto the value of zs = 3/8 given by Rankine-Hugoniot condition. There is a single jump here, so we have Fi = 1/4 and F r — 1. 
The effect of the jump condition can be interpreted geometrically as shown in the plot to the right. 



Here the superscripts +, — refer to the fact that F(z) has a positive or negative second derivative respectively. The 
superscript refers to a linear curveQ. The subscripts S mean that the solution has a shock singularity. This happens 
when F(z) becomes double valued in z and hence the solution must be resolved via a jump much like the Maxwell 
reconstruction in thermodynamics. If F starts at Fi to the left of the shock and ends at F r to the right of it the 
location of the jump is given by the Rankine-Hugoniot condition 



gCflj - G(F r ) AG 

which is simply restating the fact that the number of agents leaving the shock is equal to those that enter it. In 
geometrical terms this means that the area between the curve F(z) and z = z$ to the left of the shock is equal to 
that area to the right of the shock. 

We shall not study all the various types of solutions mentioned above. For our discussion on mergers we shall only 
need C~ and Cg . We refer the reader to FigfT] where C~ is given for P = 1 and Cg is given for T = 1. 



III. MERGER DYNAMICS 



The model presented in the previous section has a single set of rules (the set £ = {P,T,Q}) that only depends 
on the ordering of points. A straightforward generalization of this model can be to extend the rules so that we have 
various sets of probabilities and to provide a selection rule so as to define which set should be applied for a particular 
choice of agents; this selection rule should of course be a condition on the points of the agents other than their 
ordering. As we shall discuss below allowing mergers falls in a small subset of such generalizations. 



1 This case is identical to a two-agent model. Since the middle agent advances with the mean speed of the whole collection of agents, it 
is as if it is not there. 
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NO-MERGER 


LM MERGER 


LS MERGER 


MS MERGER 


MS MERGER 


MS MERGER 






L + M> S 


L + S > M 


M + S <L 


M + S> L 


M + S = L 


p 


p 2 +P9/3 


P 2 


P 2 


P 


q 


1/2 


T 


4pq/3 


pq 


q 


pq 


P 2 


p/2 


Q 


q 2 +pq/3 


q 


pq 


q 2 


pq 


q/2 



TABLE I. Winning probabilities of the three agents with points L > M > S for no-merger tournaments and for various possible 
merger cases in terms of the competitiveness parameter p. Cases where there are equalities of points can be read from this 
table using the equal likelyhood prescription described in Chapter II. For instance if M — S the probability to win for each is 
(T + Q)/2. 

A. Implementing Mergers 

The idea of mergers is only meaningful if the game is competitive, where higher points are favored. In this respect 
two agents can merge points, act as a single agent against the third one. If the mergees loose, the third player gets 
the point. If on the other hand the mergees eliminate the third player we still have to resolve a single winner. These 
considerations already hint at the necessity to construct the rules of the three-agent step in terms of a two-agent unit 
of competition. If this is so achieved, after winning against the third player the mergees can turn against each other 
and play the same two- agent unit to decide which one of them will receive the point. 

However, if the selection rule does not allow a merger, we have a generic three-agent step which we should take 
to be competitive as well to be in accord with the philosophy of the idea of mergers. Such a set of rules too can be 
established in terms of a two-agent unit of competition. In [5] it was shown that a natural way to achieve this is to 
let each three agents play a single two-agent match with each other: that is to have a tournament. All the two-agent 
matches in this tournament are decided based only on how many whole tournaments the agents have won before: 
during the tournament the tournament wins of each team (which is simply the associated points L, M and S of the 
participants) is kept constant but they accumulate match points depending on the tournament wins. The winner of 
the tournament is the agent with largest number of accumulated match points. As usual ties are decided on the basis 
of equal likelihood. 

So we have established how to implement mergers and no-mergers for each particular three-agent competition in 
terms of the rules of a two-agent unit. Having approached the problem in this way we need to use only a single 
parameter; the winning probability of the higher score agent in the two-agent unit competition. Let the rules for the 
two-agent unit be {/ > ,s} — > {p, q}: the agent with point I shall get the match point with probability p = 1 — q. This 
unit is competitive if p > 1/2. After straightforward analysis we arrive at Table|l] 

We see that various cases of mergers are possible and we need a selection rule as we have argued. The simplest and 
most natural selection mechanism which is also in accord with the overall competitive nature of our model is to let 
two agents merge if and only if both their probabilities to win the unit of competition increase with 
respect to the no-merger case. We assume that our agents are not smart and incapable of a long term strategy; 
they simply respond to an increase in the probability to win the three agent unit at hand. From the table it is evident 
that the agent with the highest point L will never find it profitable to merge neither with M nor with S. Also S 
will find it not feasible to merge with M if M + S < L. This leaves us with MS mergers with either M + S > L or 
M + S = L. Using the selection rule stated we find that if p > 3/4 mergers are feasible for those agents participating 
in i10. Analyzing further we can also discover that a merger unit is rejected if L — M > S, we remind the reader 
however that this is an interface term. 

With the analysis above we have full knowledge of the function W{x, y,y') introduced in the previous chapter. 
Other than this difference the mathematical setup of the merger model is completely analogous to the single rule 
model of the previous chapter. 

B. The Model 

The analysis of the previous section amounts to the following rules, where to be explicit at the expense of being 
redundant, all the details (including the surface terms) of the model are summarized. 



2 In fact M + S > L case only requires p > 3/5. The higher value 3/4 quoted in the text comes from the M + S = L case. This term is an 
interface contribution and thus will represent terms with higher derivatives or higher powers of the first derivative of F. Consequently 
this term is less and less important in the asymptotic future of the sytem. However in a simulation it occurs in the early stages if all 
the agents starts out with the same points so it is honest to include it in the analysis of the selection rule. 
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• L > M > S and M + S < L no merger unit with rules E = {P, T, Q}. 

• L > M > S and M + S > L merger unit with rules E' = {P', T', Q'}. 

o L > M > S and M + 5 = L merger unit with rules E" = {P", T", Q"}. 
o i = M > S no merger unit with rules {^^-, - 2^") Q}- 

o L > M = S and M + S < L no merger unit with rules jp, ^^j. 

o L > M = S and M + S > L merger unit with rules jp', T '+ Q ' , T '+ Q ' |. 

o L > M = S and M + S = L merger unit with rules jp", T "+°" , T "+°" |. 
o L = M = S equally likely unit with rules { ^ , | , ^ } . 

with the definitions P = p 2 + pq/3, T = 4pg/3 and Q = q 2 + pq/3 for no-merger units, P' — q, T' — p 2 and Q 1 = pq 
for merger units with M + S > L and P" = 1/2, T" = p/2 and Q" = q/2 for merger units with M + S = L. 

It is important to note that if only the rules S or E' are applied unconditionally during the simulation, the resulting 
distribution is in the C~ or Cg regime respectively of the single rule model. As the reader could have already guessed, 
this is somewhat evident since the set E favors the leading agent whereas the set E' is preferring the middle agent. 
Therefor when mergers are implemented the two rules are in conflict with each other. This effect is much more 
pronounced if we let the system approach the limit of full competitiveness by letting p — > 1 in which case probabilities 
converg e to E = {1,0,0} and E' = {0, 1,0} 0. We thus expect on general grounds that, if there is no prevalence of 
a single rule in the game the resulting dynamics of the system should emerge from this conflict as something that is 
neither a C~ nor a Cg solution but a new state which bears aspects of both. 

Now we have, as usual, 

r) P 

= J2J2 w ( x ~ 1 > v. Mvf' ■ ( 12 ) 

y y' 

However in contrast to the single rule model the sum above is very complicated which we do not duplicate here. In 
the hydrodynamical limit we still have, 

The particulars of it, however, are complex as expected (we have suppressed the time dependece of P's to have a 
readable expression): 

g'[F] = P' [F 2 {x) - F 2 (x/2)] + PF 2 (x/2) + 2(P - P 1 ) J ^ dy—F(x - y) 

r 2x of 

+ 2T' \F(2x) - F(x)} F(x) + 2T [1 - P(2a;)l F(x) + 2(T - T') / dy—F(y-x) (14) 

Jx dy 

f°° dF 

+ Q'[1 + F 2 {x) - 2F{x)F{2x)] - 2Q [1 - F(2x)} F(x) + 2(Q - Q') J dy—F(y - x) 

As a check we see that letting P' = P and T" = T we recover the single rule model. This equation has all the 
complicating adjectives one can attach, the most important being non-locality, and a direct approach as that of 
solving for characteristics is not obvious. Nevertheless the dependence on x/2 and 2x in (fl"4|) are suggestive and have 
a rather simple interpreation; an agent with point x is protected against mergers of two others if the points of 



3 One may object to this limit by observing that the probability to win for the agent with the lowest point remains zero and thus this 
player has no incentive to participate in a merger with the middle agent. However letting p = 1 — 9 and expanding the probabilities to 
first order we see that Q = 0/3 and Q 1 = 9. Thus the lowest lagger triples its chances to win in merging no matter how close to zero its 
winning chances are. 
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FIG. 2. Results of the numerical simulation for various values of the competitiveness parameter p. For a detailed analysis we 
refer the reader to the text. 



those are both smaller than x/2 and similarly an agent with point 2x is protected against mergers of two agents with 
points less than x. 

One can also show that the ansatz F(x, t) = F(z = x/t) of the single rule model is still applicable her43. The only 
concern could be the integrals but they are easily transformed accordingly. For instance, 

j dy—F{x -V)= J dtF'(()F(z - , (15) 

where £ = y/t. 



C. Numerical Analysis 



In view of the obvious difficulty of the equation we resorted to numerical simulation with the hope that it will give 
us clues about the system. We have performed extensive numerical studies on the system. To be spesific we simulated 
the microscopic system from p = 0.750 to p = 1.000 in steps of 0.001. For each of these the run was consisted of 
N = 10 6 agents and we ended the simulation when on average each agent had played about 6.5 x 10 6 games meaning, 
in our normalization, that we stopped the run when time variable t is 2 x 10 . For all these cases the initial condition 
was F(x,0) = Q{x). 

To start the exposition of numerical results we invite the reader to analyze FigJ5] One important aspect one can 
immediately realize is that upto about p = 0.76 the resulting distribution is exactly the same as a pure C~g game. This 
can only happen if the game is dominated with mergers and no-mergers never happen; so essentially we have a single 
rule game. Let us call the location of the jump zr and denote point where F vanishes as zl- Mergers can gloablly 
dominate if and only if zr > 2zl which simply implies that for any choice of three agents M + S > L. So for this case 
we can use the results of the single rule model to predict when zr = 2zl which yields po — (V13 + l)/6 « 0.76759 



4 This is a consequence of the fact that we have W(ax, ay, ay') = W(x, y, y') for either the single rule or the merger model. This in turn 
means that the equation has x ax and t — > at symmetry. 
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FIG. 3. The qualitatively stable (see text) solutions of the merger model for various number of bunches. 



in accordance with the simulations. After this point we see that the number of agents in the vicinity of zl increase 
0. This happens because after the transition the agents accumulated at the jump discontinuity becomes protected 
againts mergers of the agents around Zl and some agents slide down the slope. These agents able to merge no more 
against the the bunch at zp lag faster than before. Increasing p emphasizes this effect. Later on during the excursion 
to higher values of p, to the right of Zl a flat plateau appears. Let us call the position of the left tip of this plateau as 
zp. The emergence of this plateau coincides with the condition zp — 2zp at around p = 0.801. This means that the 
bunch condensed on the discontinuity at zp becomes completely protected from mergers to the left of Zp. Furhtermore 
since zp > 2zp the games played among agents to the left of zp are all merger dominated and thus locall^ of type 
Cg. Increasing p further we see that another condensation of agents occurs at zp, which itself shifts to the left as p 
takes on higher values. We have turned around full circle and the solution consists of a right bunch condensed at zp 
and a left bunch sitting at values less than or equal to zp. The local structure of this left bunch is similar to the the 
whole structure when p was less than pq. As the reader could have guessed a further excursion to higher values of p 
repeats the structure and we end with a seemingly self similar pattern. Our numerical analysis indicates that as p 
approaches this pattern repeats itself ad infinitum. 

To recapitulate we refer to Fig(3] where we have provided representative solutions relevant to our discussion. There 
the solutions we have presented are those that have a qualitative stability, that is the location and the heights of the 
shocks and the form of the leftmost group will alter as we change p but the number of shocks will be constant for a 
while during such an excursion. We shall call these shocks bunches and label them as B n with n = 0, 1, 2 • • • were 
n = representing the rightmost bunch. The leftmost bunch is not a pure shock but a combination of a shock and a 
rarefaction wave since it is locally a C<t solution. The transitions between those solutions are complicated and bear 
the full complexity of the differential equation which is somewhat unpenetrable. As can inferred from Fig. [31 the most 
important aspects of these solutions are 

• The games among players of a single bunch B n are all mergers, since the points of these agents are localized 
around a particular value. 



• All players in a bunch B n are protected from mergers of two players in Ufe>n 
each bunch is larger than twice that of those to the left of it . 



-Bfc, since the localized points of 



5 Needless to say zr and zl are changing as p increases. As can be expected on general grounds zr moves to the right and zl moves to 
the left. 

6 That is if we consider the games only between agents in the left bunch. 
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D. Self similar behaviour of the model 



In view of the discussion above, for values of p larger than about 0.81, we can separate the equation into two 
parts, one valid near the condensation of the rightmost bunch (the forrunncr agents) and the rest. All we need is the 
interaction of these two parts. From the numerical study we know that the rightmost bunch is completely protected 
from mergers of two agents not in this bunch and that the games within this bunch are merger dominated. Let us 
start with the generic form of the eqution 



As we have mentioned for the merger model we can recast this in terms of z = x/t 



— (- z + g'[F}) = (17a) 

g'[F] = J J dydy'W(x,y,y')^^ , (17b) 

where the integrals are over the whole domain of points. 

Let us denote the rightmost bunch as B and all the agents not in this bunch as B^. The points of the agents in 
B are coalesced near a value Xo(t) which is at least twice as large as the largest point of the agents in B^. Therefore 
the equation for B can be written as 



dF 0F Q f f dFdF 
where in evaluating W we should remember that x is near xq. We end up in 



dF 
dt 



dF 
dx 



(19) 



with 



Q'o =P'(F - $x) 2 + 2T'(F - $0(1 - F ) + Q'(l - F ) 2 (20) 
+2P'$ 1 (f 1 - $i) + 2T'$i(l - F ) 
+P$ 2 

The value $i represents the total number of agents in B^. That is, Fq starts at $i to the left of Bo and ends at 1 
to the right. The first line above are the self games of the rightmost bunch. The second line represents games where 
two agents are in B and one in B^. The third line represents games where only one agent is selected from B a . 
For values of x < x /2 we are in B^. Denoting the cumulative function in this region as FL(x,t) we get 



with 



dF L 
dt 



dF L 

dx 



(21) 



g' L =J J dydy'W(x,y,y>) 



dF L 0F L 



dy dy' 

+2T(1 - *i)F L + 2Q(1 - $i)(*i - F L ) 
+Q'(l-$i) 2 



(22) 



(23) 



Now Fl starts from and ends at $i and the second (third) lines in the equation above represents picking two (one) 
agents from Bl- The first line represents the self games of all the agents that are in Bl- 
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There is resemblance to self-similarity in (|2"2"|) ; the first line looks like the original equation but the interaction terms 
with the bunch Bq spoils this as Bl can still gain points via these terms. However in the extreme competitiveness 
limit p —> 1 these terms are absent since T, Q and Q' all vanish. That is, agents in Bl will only gain points againts 
themselves and will simply remain idle during any competition with the bunch Bq. Conversely bunch Bq will use Bl 
as a definite source of points. In this limit (j2"2")l becomes 

Q' L = j J dydy'W(x,y,y')^^ (24) 

Let us recall however that the maximum value Fl can take is $i. Defining Fl = Fl/^i and using the ansatz z = x/t 
for a solution the equation becomes 



dF L {z) 
dz 



(25) 



This has exactly the same form as the original equation if we let z — > <&\z which one could interpret as scaling of 
x. However unless the initial condition can be partitioned this way we can not say that the solution will resolve itself 
into a self-similar shape. The inital data we use F(x, t) = 0(x) can be partitioned this way because Q(ax) — Q(x). In 
conclusion if F(x,0) — Q(x) we expect self-similar behaviour in the solution when p = 1 via the scaling we 
have described above. The procedure of extracting the rightmost bunch is somewhat similar to renormalization 
procedure and the decoupling mechanism in field theory where after integrating out high energy degrees of freedom 
we end up with a new theory. If the original theory is said to be non-renormalizable the new theory is different. If 
otherwise the new theory is similar in form to the original except quantities in it like fields, coupling constants etc. are 
scaled it is called a renormalizable theory. We see an analogy here; the extraction of the rightmost bunch yields the 
same form of equations for p = 1 and different otherwise. Thus within this sense we can say that for p = 1 and with 
F(x, 0) = Q(x) the model is renormalizable. The effect of this renormalization yields the scaling of F via F — > $iF 
and that of z via z — > z§\. 

Given these circumstances we can repeat the same structure infinitely many times where we end up with infinitely 
many bunches B n localized around z n containg f n = $ n +i — $„ agents. And as expected at each step the scaling 
should be via the same number $i = ip if there is to be self-similarity at all. The protection of B n againts mergers 
of any two players in -Bfc> ra mandates that z n > 2z n +\. Reiterating this procedure we find the equation obeyed for 
agents in B n to be 

dF 

[-z - 2F n (z) 2 + 2F n (z)$ n + $ 2 n+1 ] = . (26) 
As expected this can only be resolved via a shock, the location of which is found via the Rankine-Hugoniot condition 

*n = ^K+i + <MW + (27) 
Since at each step we scale with the same number <p we have 

$n = ■ (28) 

Which in concert with the protection mechanism mentioned means that <p < l/v2 and implies the following 

z n = z oV 2n , (29a) 
zo = + <P + 1) • (29b) 

We are one equation away from a solution. A concept we may use is the stability of the solution at large times. 
For instance let us focus on the second bunch B\ . To the right of this bunch there is Bq containg 1 — ip agents and to 
the left of it there are infinitely many bunches accomodating (p 2 elements. Now if the solution is stable the location 
of Bi in terms of the variable z does not change in time. That is the points they loose against Bq is the same as the 
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FIG. 4. The solution to the merger model in comparison to that of the single rule model both for p = 1. The mean speed of 
agents is shown with a vertical line. 



points they gain from all the agents to the left. This can only happen if these two regions have the same number of 
agents yielding tp 2 — 1 — ip. One can similaryly argue as follows; how can the bunch Bq knows that it is the leading 
bunch in a self-similar pattern? Let us assume the existance of a further bunch B—\. Using the scaling we expect 
the number of agents in this bunch to be 1/cp — 1. Now the location of Bq will be stable if and only if the games lost 
to B-i only by Bq equals the games won against agents in Uk>oBk again only by Bq. This can only happen if the 
number of agents in _B_i equals the number of agents below Bq meaning l/(p — 1 = (p. These two considerations are 
equivalent and allows us to find 



V5-1 

ip = — - — , (30a) 



zo = I ■ (30b) 



Thus (p is the reciprocal of the Golden Ratio and zq is just twice the value of the mean speed of the entire system. So 
a player in the rightmost bunch is, in the mean, on the verge of being protected from mergers of any two randomly 
picked agents from the entire collection. The simulation results for which F(x, 0) = Q(x) are in very good aggreement 
with this theoretical prediction. 

It is interesting to contrast the model with mergers and the single rule model without mergers. The comparison is 
in Fig. 21 One important aspect we realize is that the effect of implementing mergers does not in general affect the 
global behaviour of the problem; the model with mergers is like a discretization of the curve F(z) = y/z so overall 
competitiveness is still there. However the local behaviour is completely different; we have stratification of agents. 
It is rather peculiar to observe this behaviour when we, in effect, increased the overall competitiveness of the model. 
That is, one can say that the environment is now subject to more coercive rules of competition which in contrast 
leads to condensation of agents around particular values of point gain. As we have mentioned before this effect is a 
consequence of the conflict between mergers living in the Cg regime and no-mergers being in the C~ class of solutions 
were the simulations run unconditionally as a merger or no-merger single rule model respectively. 



13 



E. Restricted Mergers 



The merger model we have presented represents a very coercive environment of competition and it is not readily 
susceptible to analytical study for arbitrary p. The main reason for this is its high non-locality. This non-locality is 
present because the two lowest laggers are allowed to merge in all cases even when M + S = L + 1. One could wish to 
contemplate other schemes of mergers where this effect is less pronounced. One way to do this is to regulate mergers. 
Here we present a model which is the most restricted. Let us remember that at time t the theoretical maximum point 
is just t. Now let us pick three agents and order their points as L > M > S and let us allow mergers only if S > t/2; 
this makes sure that the competition is a merger unit since we always have L < t and thus it is always true that 
M + S > 2L. In such an approach mergers will be represented in W via a term like 0(x — t/2) wich will become 
0(z — 1/2) for the asymptotic behaviour where as usual z = x/t. The equation becomes 



dF 



-z + G'(F)} = 



(31) 



with 



Q'(p) = J PF2 + 2TF(1 + 0(1 - F f for z < I/ 2 f32 ) 

1 ' \PF 2 + {P' -P)(F-F) 2 + 2TF{l-F) + 2T'(l-F)(F-F) + Q / {l-F) 2 for z > 1/2 

and with F = F{\/2). This equation is local and thus can be studied analytically in much the same way as the single 
rule model. Here we only present the solution for p = 1 to compare it with the unrestricted merger model. 




(33) 



where we have 



_ 2a/2- 1 

Z r - ^ ' 

Which is in very good agreement with numerical simulations. As we see the solution for z < 1/2 is the same as that 
of the single rule model. The comparison for the full range of z is given Fig. [5j As the reader could have guessed the 
stratification effect is non-existent but the tendency of the curve to go to that regime had mergers were unrestricted 
is apparent. 



IV. DIGRESSION ON INITIAL CONDITIONS 



Let us remember that the generic form of the equation governing the dynamics of three agents games is 



dF dF 



For the single rule model of Chapter I. the integrals resolve into a simple polynomial of F. The reason for such a 
simplification is, for the model mentioned, the fact that W being only a function of the ordering of points containing 
only terms like &(L — M)Q(M — S) and thus obeying 



W(ax, ay, ay') = W(x, y, y') 
W(x -b,y-b,y' -b) = W(x, y, y) 



(34a) 
(34b) 
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FIG. 5. The solution to the restricted model case in comparison with the curve for single rule model both for p = 1 

The first of these equations means that the equation will be invariant under the combined transformations x — > ax 
and t — > at which makes it possible to assume an ansatz F{x/t) since F(x, 0) = Q(x) also obeys this symmetry. 

Now let us shift F by a constant </> such that F = F + (f> with <fi = — (T — Q)/(l — 3T). Since the single rule model 
obeys (|34bj) one can also perform a Gallilean transformation^] on the independent variables of the form t = t and 
x = x — zot. Choosing r = Q — (T — Q) 2 / (1 — 3T) the equation becomes 



which can be recast as 



G(F) = (1 3r) F 3 . (35b) 

where G(F) is strictly concave for (1 — 3t) > 0. If on the other hand (1 — 3i) < one can make G(F) strictly concave 
by doing x —> —x. 

After these transformations the initial condition is changed into F(x, 0) = Q(x) + 4> which still represents a Riemann 
problem. Now it is a known fact that the solutions to equations of type (|35p converge in the infinite time limit to the 
solutions of the Riemann problem^ if the initial condition obeys F(x,0) = Fl for some x < xl and F(x,0) = for 
some x > xr where Fl and Fr are constants. Now the dependent variable F is the cumulative of a globally conserved 
quantity; the number of agents. Therefor for a generic initial distribution of points we have F(x, 0) = for x < 
and F(x, 0) = 1 for some x > xr. We thus infer that for the single rule model the time asymptotics of F is 



7 Under the Gallilean transformatio alone the equation is transformed into [—2 + zq + G'(F)\. That is the z variable which in principle 
represents the speed of agents is shifted by a constant. 

8 See for instance [|, [j| and p]} 
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independent of the initial conditions^]. These considerations also apply to the restricted merger model we have 
studied since in principle it has the same general form as the single rule model. 

For the merger model without restrictions we still have the symmetry in (|34al) which allows us to make the F(x/t) 
ansatz if the initial condition obeys F(ax, 0) — F(x,0). Unfortunately the shift symmetry in (|34b|) is absent because 
mergers are implemented via terms of the type 0(L — M )<d(M — S)0(M + S — L). This makes the equation highly non- 
local and in particular the Gallilean transformations will take it to an entirely different form. None of the theorems 
presented in the mentioned papers above hold and one would expect a strong dependence of the time asymptotics on 
the initial data; an observation which we have substantiated with numerical simulations. 



V. DISCUSSION 

The unrestricted merger model we have presented has interesting properties. The most important being the stratifi- 
cation of the entire society of agents. The bunch that has the largest number of agents (this number is (1 — tp) w %32) 
is also the bunch with the largest rate of point gain. However this point gain is only two thirds of the maximum 
possible rate. On the other hand all the agents other than the first bunch are earning slower than the mean rate. The 
number of agents living below this mean is slightly lower in the single rule model where mergers are not implemented. 
Furthermore from our solution it is clear that the number of agents in a bunch f n satisfies /„ = f n —i — fn+U that is 
the number of agents in a bunch is like a derivative in the sense of the bunches. 

It is tempting to speculate that the merger model we have presented could have applications to natural or social 
phenomena. The stratification phenomenon being present in various systems. For instance one could argue that a 
bunch, in essence, represents a single entity, the number of agents in it simply meaning that it has more activity 
in taking part in games. From this perspective one may interpret the merger model as one of explaining monopoly 
formation after a period of competition between companies. Stratification is also present in natural systems. Another 
tempting interpretation could be the stratification of the collection of entire living species in terms of their genetic 
material. If there is a competition mechanism complexified with mergers like the one described in this work one could 
hope to gain qualitative understanding of the formation of different strata of living organisms. Agents could be units 
of genetic material and the competition could be for taking part in the genetics of a (new) species. 

The emergence of the stratification mechanism can be interpreted in the following way. The rules of mergers yield 
a complicated and non-local system. The equations are so complicated that they can not be resolved via smooth 
functions and the only possible escape is the formation of various shocks; there must be a solution since we are 
simply simulating a Monte-Carlo system with a well defined microscopic model. We believe this to be true for other 
systems of conservation laws, coming from well defined microscopies, where the equations become non-local. 
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